A log-Birnbaum-Saunders Regression Model with 

Asymmetric Errors 



o 

(N 

Oh: 

in 



C/3 



> 

oo 
oo 

00 

p 

o 
o 



X 



Artur J. Lemonte 

Departamento de Estatistica, Universidade de Sao Paulo, Sao Paulo/SP, 05508-090, Brazil 



Abstract 



The paper by 



Leiyaetal 



(120101 ') introduced a skewed version of the sinh-normal distribu- 
tion, discussed some of its properties and characterized an extension of the Birnbaum-Saunders 
distribution associated with this distribution. In this paper, we introduce a skewed log-Birnbaum- 
Saunders regression model based on the skewed sinh-normal distribution. Some influence meth- 
ods, such as the local influence and generalized leverage are presented. Additionally, we derived 
the normal curvatures of local influence under some perturbation schemes. An empirical applica- 
tion to a real data set is presented in order to illustrate the usefulness of the proposed model. 

Key words: Birnbaum-Saunders distribution; fatigue life distribution; influence diagnostic; max- 
imum hkelihood estimators; sinh-normal distribution; skew-normal distribution. 



1 Introduction 

The two-parameter Birnbaum-Saunders (BS) distrib ution, also known as the fatigue life distribution, 



was introduced by 



Bimbaum and Saunders 



(1969a 



H). It 



was originally derived from a model for a 



physical fati gue process wher e dominant crack growth causes failure. A more general derivation was 
provided bvlDesmondl (119851) bas e d on a bi ological rnodel a nd relaxing several of the assumptions 
made by iBimbaum and SaundersI (Il969a|) . lOesmondl (Il986|) investigated the relationship between 
the BS distribution and the inverse Gaussian distribution. The author established that the BS distri- 
bution can be written as a mixture equally weighted from an inverse Gaussian distribution and its 
complementary reciprocal. 

The random variable T is said to have a BS distribution with parameters a,r] > 0, say BS(a, r]), if 
its cumulative distribution function (cdf) is given by F{t) = $(f ), t > 0, where $(■) is the standard 
normal distribution function, v = p{t/r])/a, p{z) = z^l"^ — z^^/"^ and a and rj are shape and scale 
parameters, respectively. Also, rj is the median of the distribution: F{ri) = $(0) = 1/2. For any 
constant > 0, it follows that kT ~ BS(a;, ky) . It is noteworth y that the reciprocal property holds 



for the BS distribution: T ^ ~ BS(a;,?7 ^); see ISaundersI (|1974|) . The BS distribution has received 
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considerable attention over the last few years. 



KunduetaL 



(|2008|) discussed the shape of the hazard 



function of t he BS distribut i on. R esults on improved s t atistical infe rence for the BS distribution are 



discussed in 



Wu and WongI (120041) and 



Lemonte et al. 



(2007, 



tensions of the BS d istrib ution are presented i n iDfaz-Garcia and Leival (|2005b . 



Guiraud et al. 



(I2OO9I) and 



Castillo et al. 



2008h. Some g e neralizations and ex - 



Gomes et al. 



ies (see, fo r exam ple, iBalakrishnan et al.L l2007h and outside t his field; see 



Leiva et al. 



(l2009h. T his distribution has been ap plied in reli a bility stud- 

( 2008 ) and 



Leiva et al. 



(I2009D . Additionally, based on the BS distributio n, iBhattil (120101) introduced the BS au- 
toregressive conditional duration model. IXu and TangI (120 10|) presented estimators for the unknown 
paramete rs of the BS d istribution using reference prior. 



FromRiec 



19891), if 



Z = v -\ sinhf —] 

a \ a J 



N(0,1), 



(1) 



then Y has a four-parameter sinh-normal (SHN) distribution, denoted by F ~ SHN(a, 7, a, u), where 
z/ G 9ft and a > are the shape parameters, and 7 £ 'St a nd a > correspond to the location and 
scale parameters, respectively. According to iRieckl (Il989() . the parameter u is also the noncentralty 
parameter. If z/ = 0, the notation is reduced simply by F ~ SHN(a, 7, a), and this distribution has 
a number of interesting properties. For example, it is symmetric around the mean E(F) = 7, it is 
unimodal for a < 2 and bimodal for a > 2 and if ~ SHN(a, 7, a), then = 2{Ya — 7)/ (aa) 
converges in distribution to the standard normal distribution when a — 0. If F ~ SHN(a, 7, a = 2), 



then T = exp(F) follows the BS distribution with shape parameter a and sea 



Leiva et al. 



e para meter rj = 
(l2010l) . the SHN 



exp(7), i.e. T = exp(F) ~ BS{a,ri). For this reason, according to 
distribution is also called the log-Birnbaum-Saunders (log-BS) distribution. Additionally, according 
to these authors, the SHN and BS models correspond ing to a logarithmic distribution and its associated 
dis tribution, respectively ( Marshall and Olkin . 2007[ Ch. 12). 



Rieck and Nedelman 



^1991|) introduced a log-BS regression model based on the SHN(a;,7,2) 
dist ribution. Their r e gression model has been studied by sev e ral authors. Some impo r tant references 



(2010 


), 


Xiao et al. 


(2010 


) and 


Cancho et al.l ( 


2OI0I). 



Desmond et al. 



(l2007h . 

among others. General 

reg ression model introduced bvlRieck and NedelmanI (| 199 11) are presented in 



Lemonte et al. 



(l2008l) . 
izations of the log -BS 



Xi and Wei (12007 



and 



4) 

Lemonte and Cordeirol (|2009h . 
Leiva et al.l ( 201o|) introduced a skewed SHN distribution by replacing the standard normal distri- 



bution in equation ([T]) by the skew-normal (SN) distribution (|Azzalinel Il985|) . i.e. they consider the 
random variable 



— sinh ( — ) 

a \ a J 



SN(A), 



where A G 9ft is the shape parameter which determines the skewness. Now, the notation used is F ~ 
SSN(a, 7, a, z/, A). From now on, we shall consider u = and a = 2 and hence the notation is given 
by F ~ SSN(a, 7, 2, A). The random variable T = exp(F) follows the extended Birnbaum-Saundres 
(EBS) distribution, with shape parameters a > and A G 9ft, and scale parameter i] = exp(7). Now, 
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the notation is T = exp(y) ~ EBS(a, r], A). 

Let T ~ EBS(a, r], A). The density function ofY = log(r) is given by JLeiva etallboiol) 



cosh 



(L 



7 



a 



\ 2 



sinh 



(y. 



7 



a 



\ 2 



$ ( — sinh 
\ a 



y-1 



where 0(-) is the standard normal density function, and, as before, we write Y ~ SSN(a, 7, 2, A). 
The sth (s = 1,2,...) moment of Y can be written as 



E(n = 2'=2Z7'-'cfc(«,A), Cfe(«,A) 

fc=0 

Thus, the mean of Y is given by E(F) = 7 + c{a, A), with 



/oo 
{sinh~^(aw/2)}V(w)$(Aw)c/w. 
-00 



c(a. A) 



/oo 
{sinh"^(att;/2)}0(w;)$(Aw)rfM;. 
'OO 



Plots of the SSN(q;, 7, 2, A) distribution are illustrated in Figure[T]for selected parameter values. 
The chief goal of this paper is to introduce a skewed 



SSN(a, 7, 2, A) distribution, recently proposed by 



LeivaetaL 



og-BS regression model based on the 
(120101) . The proposed regression model 
is convenient for modeling asymmetric da ta, and it is an alternative to the log-BS regression model 
introduced by iRieck and NedelmanI (Il99l|) when the data present skewness. The article is organized 
as follows. Section [2] introduces the class of skewed log-BS regression models. The score functions 
and observed information matrix are given. Section [3] deals with some basic calculations related with 
local influence. Derivations of the normal curvature under different perturbation schemes are pre- 
sented in SectionlH Generalized leverage is derived in Section[5l Section [6] contains an application to 
a real data set of the proposed regression model. Finally, concluding remarks are offered in Section|71 



2 Model specification 

The skewed log-BS regression model is defined by 

yi = xj(3 + ei, z = l...,n, (2) 

where yi is the logarithm of the ith observed lifetime, Xi = (xn, . . . , Xip)~^ is a vector of known 
explanatory variables associated with the ith observable response yi, f3 = . . . , (3pY is a vector 
of unknown parameters, and the random errors Ei ~ SSN(q;, —c{a, A), 2, A) that corresponds to the 
regression model where the error distribution has mean zero. Thus, we have y,; ~ SSN(a;, xjf3 — 
c(a. A), 2, A), with E(?/j) = xj^, for z = 1, . . . , n. 

The log-likelihood function for the vector parameter Q = a, A)^ from a random sample 
y = (yi, . . . , yn)~^ obtained from Q can be expressed as 

n 

i{0) = Y,UO), (3) 

i=l 



3 




a = 1 .5 and Y = a = 2 and y= 




a = 2.5 and y=0 a = 3.5 and y=0 




Figure 1: Plots of the density function of the SSN distribution for some parameter values. 
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where £,(0) = - log(27r)/2 + hg{^a) - ^1/2 + log{$(Ae.2)}, 

in = 6i (^) = - cosh , ^i2 = 62 (^) = - smh — 

a \ 2 / a \ 



- xjl3 + c{a, A) 



for i = 1, . . . ,n. The function i(0) is assumed to be regular dCox and Hinkleyl 1 19741 Ch. 9) with 
respect to all /3, a and A derivatives up to second order. Further, the n x p matrix X = (xi, . . . , ajn)^ 
is assumed to be of full rank, i.e., rank(X) = p. 

By taking the partial derivatives of the log-likelihood function with respect to (3, a and A, we 
obtain the components of the score vector Ue = {Uj ,Ua,Ux)^ ■ We have Uf3 = X^s, where 
s = {si,...,Sny withsj 



= - UlU - Aea0(A6i)/$(Aei2)}/2, 

0(A^i2) { Cain 



ii2 

a 



n 
i=l 



1 ^{K^2 

^1 _^ 1 ^ 0(A62 



where 



Ca = Ca(a, A) 

Ca = ca(«, A) 



4 



■{Xcxiii + 2ei2}, 



w(4 + a^w'^Y^/'^(t){w)^{\w)dw, 
w smh.~^ {aw /2)(f){w)(^{Xw)dw . 



Setting these equations to zero, Ue = 0, and solving them simultaneously yields the MLE = 
(/3^, a, A)^ of = (/3^, a, A)^. These equations cannot be solved analyticall y and statistical soft- 
ware can be used to solve them numerically. For example, the BFGS method (see. lNocedal and Wright! 



1999; 



Press et al. 



20071) with analytical derivatives can be used for maximizing the log-likelihood 



function i{6). Starting values /3^'^\ a'-^^ and A*^°-' are required. Our suggestion is to use as an 
initial point estimate for /3 the ordinary least squares estimate of this parameter vector, that is, 
13 = {X^ X)^^ X^y. The initial guess for a we suggest is where 



a 



— sinh^ 
n ^ 



2fyi- xjf3 



i=l 



We suggest A*^*^-* = 0. These initial guesses worked well in the application described in Section [6l 

The asymptotic inference for the parameter vector 6 = {f3^ , a, A)^ can be based on the nor- 
mal approximation of t he MLE of 6, = (/3^,a, A)^. Under some regular conditions stated in 
ox and Hinkleyl (|l974L Ch. 9) that are fulfilled for the parameters in the interior of the parameter 
space, we have ~ Mp+2{&, S^), for n large, where ~ means approximately distributed and Eg 
is the asymptotic variance-covariance matrix for 6. The asymptotic behavior remains valid if is 
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approximated by —Lgg, where —Lgg is the {p + 2) x {p + 2) observed information matrix evaluated 



at d, obtained from 



where 





Li3l3 


Lfia 


Lf3\ 




-X^VX -X^h 


-X^b 


Lee — 


Laf3 


Laa 


LaX 




-h^X tr(Ki) 


tr(K2) 




Lxp 


L\a 


Lxx_ 






tr(K3)_ 


diag{wi, . . . , 






= diag{fc 




= diag{/c 


K3 = diag{A;i3, . . 






h = 


(/ll,...,/l„)^, b 





j2; 



All the quantities necessary to obtain the observed information matrix are given in the Appendix. 



3 Local influence 



The local influence method is recommended when the concern is related to investigate the model 
sensibility under some minor perturbations in the model (or data). Let lj G f2 be a /c-dimensional 
vector of perturbations, where O C Sft'^ is an open set. The perturbed log-likelihood function is 
denoted by i(0\uj). The vector of no perturbation is cjq E fl, such that i(0\ujQ) = i{0). The 
influence of minor perturbations on the maximum likelihood estimate 6 can be assessed by using the 
log-likelihood displacement LD^^ = 2{£{0) — 1(0,^)}, where 0,^ denotes the maximum likelihood 
estimate under i(0\uj). 

The Cook's idea for assessing local influence is essentially to analyse the local behavior of LD^^ 
around by evaluating the curvature of the plot of LD^^^^ad against a, where a G 3? and d is a 
unit norm direction. One of the measures of particular interest is the direction d^^x corresponding 
to the largest curvature Cd^^^. The index plot of dmax may evidence those observations that have 
considerable influence on LD^^ under minor perturba tions. Also, plots of dmax against covariate 
values may be helpful for identifying atypical patterns. ICookl (|l986h shows that the normal curvature 
at the direction d is given by 

Cd{e) = 2\d''A'^L,^Ad\, 

where A = d'^i{6\(jj)/d0duj^ and —Lee is the observed information matrix, both A and Lee are 
evaluated at 6 and cvq. Hence, is the largest eigenvalue of = —A^LggA and dmax is 

the corresponding unit norm eigenvector. The index plot of cZmax for the matrix B may show how to 
perturb the model (or data) to obtain large changes in the estimate of 6. 



Assume that the parameter vector is partitioned as 6 
02 are pi and p — pi, respectively. Let 



[6j ,6JY . The dimensions of 6i and 



Le, 01 



'0102 



Lei02 

-C'02 02 



6 



where Lg^e, = dH{e)/deideJ , Lg^g^ = dH{e)/ddideJ and Lg^g^ = d'^ i{e)/ 862861 . If the 
interest lies on 61, the normal curvature in the direction of the vector d is Cd-^g^ (6) = 2\dJ A.^ {Lgg — 
Adl, where 



■122 



O2 02 



and dm ay;fli here is t he eigenvector corresponding to the largest eigenvalue of Bi - 
1,22)^ (ICookl . ll986h . The index plot of the d^^^^-g^ may reveal those influential elements on 61. 



4 Curvature calculations 



Next, we derive for three perturbation schemes the matrix 

. 8''m\u)) 



868uj' 



A/3 

A, 
Aa 



considering the model defined in Q and its log-likelihood function given by ([3]). The quantities 
distinguished by the addition of are evaluated did = (/3^, S, A)^. 



4.1 Case- weights perturbation 

The perturbation of cases is done by defining some weights for each observation in the log-likelihood 
function as follows: 

n 
i=l 

where uj = (wi, . . . , a;„)^ is the total vector of weights and wq = (1, . . . , 1)^ is the vector of no 
perturbations. After some algebra, we have 

A/3 = X'^S, A„ = (oi, . . . ,a„), Aa = (ci, . . . ,c„), 

where S = diag{si, . . . , 

c. = --[Ui.2 - ^1 + ^^p;^(262 + Aca6i), 

fori = 1, . . . ,n. 
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4.2 Response perturbation 

We shall consider here that each yi is perturbed as = Vi + cOiSy, where Sy is a scale factor that 
may be estimated by the standard deviation ofy. In this case, the perturbed log-likelihood function is 
given by 

n 1 " 

i{e\u^) = -^log(87r) + J]log(ea.J - ^J^^L.^ 

i=l i=l 

where = ^ilw^i0) = 2a~^ cosh{[y i,^ - xj f3 + c{a, X)]/ 2), ^i2w^ = iiiwA^) = 2«"^ sinh([?/i^ - 
xjf3 + c{a, A)]/2) and ojo = (0, . . . , 0)^ is the vector of no perturbations. Here, 

= SyX^V, A„ = Syh'^, Aa = Syb^. 



4.3 Explanatory variable perturbation 

Consider now an additive perturbation on a particular continuous explanatory variable, namely Xj, by 
making Xijw = Xij + WjS^, where is a scale factor that may be estimated by the standard deviation 
of Xj. This perturbation scheme leads to the following expression for the log-likelihood function: 

n 1 " 

i{6\u) = log(8vr) + J]log(e.i.J - 2 E^'^.., 

i=l i=l 

where ^n^,^ = ^ii^^{0) = 2a-^ cosh{[yi - xj^f3 + c(a, A)]/2), ^i2«;2 = 62u,2(^) = 2a"^ sinh([?/i - 
xl,^ + c{a, A)]/2), with Xi^ = {xn, Xij^, Xip)^ . Here, cjq = (0, . . . , O)'^ is the vector of no 
perturbations. Under this perturbation scheme, we have 

A/3 = -s^^jX^V + s^CjS^, A„ = -sJ3fi7, Aa = -sji^i^, 

where Cj denotes a p x 1 vector with 1 at the jth position and zero elsewhere and 13 j denotes the jth 
element of /3, for j = 1, . . . , p. 



5 Generalized leverage 



Wei et al. 



(119981) . which is defined 



In what follows we shall use the generalized leverage proposed by[ 
as GL{6) = dy/dy^, where 6 is an s-vector such that E{y) = ^i{6) and 6 is an estimator of 0, with 
y = n{6). Here, the (i, /) element of GL{0), i.e. the generalized leverage of the estimator 6 at (i, I), 
is the instantaneous rate of change in ith predicted value with respect to the Ith response value. As 
noted by the authors, the generalized leverage is invari ant under reparameterization and observations 
with large GLij are leverage points. Iwei et aL ( 1998 ) have shown that the generalized leverage is 
obtained by evaluating 

GL(e) = Df,(-Lf,f,)-'L, 



Dg{-Lg0) 



ate = 0, where Dg = dfi/dO^ and Lgy = dH{9)/d0dy 
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After some algebra, we have that 



X 



and 



Thus, from these quantities, we can obtain the generalized leverage. 



6 Application 



In this section we shall illustrate the usefulness of the proposed regression model. The fatigue pro- 
cesses are by excellence ideally modeled by the Birnbaum-Sau nders distribution d ue to its genesis. 



Chan et al. 



( 20081) . These data con- 



We consider the data set given in lMcCooll(|1980|) and reported in 
sist of times to failure (T) in rolling contact fatigue of ten hardened steel specimens tested at each of 
four values of four contact stress (x). The data were obtained using a 4-ball rollin g contact test rig a t 
the Princeton Laboratories of Mobil Research and Development Co. Similarly to 
we consider the following regression model: 



Chan et al. 



mm. 



Ui = 1^1 + h ^og{xi) + Ei 



1, 



40, 



where yi = log(Tj) and £i ~ SSN(a, — c(a, A), 2, A), for i = 1, . . . , 40. All the computations were 



done using the Ox matrix programming language (lDoomikLl2006|) . Ox is freely distributed for aca- 
demic purposes and available at http://www.doornik.com. 



Table 1: Maximum likelihood estimates. 



log-BS skewed log-BS 

Parameter Estimate SE Estimate SE 

Ji 0.0978 0.1707 0.1657 0.1759 

^2 -14.1164 1.5714 -13.8710 1.5887 

a 1.2791 0.1438 2.0119 0.3487 

A — — 1.6423 0.5679 

Log-likelihood -61.62 -58.68 

AIC 129.24 125.36 

BIC 134.31 132.12 

HQIC 131.07 127.80 



Table [H lists the MLEs of the model parameters, asymptotic standard errors (SE), the values of 
the log-likelihood functions and the statistics AIC (Akaike Information Criterion), BIC (Bayesian 
Information Criterion) and HQIC (Hannan-Quinn Information Criterion) for the skewed log-BS and 
log-BS regression models. The SE of the estimates for the skewed log-BS model were obtained using 
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(a) 



(c) 



(b) 



(d) 



Figure 2: Index plots of |<irnax| for under case weighting (a), response (b) and covariate (c) pertur- 
bations, and generalized leverage (d). 

the observed information matrix given in Section[2l while the SE of the estim ates for the log-BS m odel 



were obtained using the observed information matrix given, for example, in lGalea et al.l (|2004|) . The 
estimatives of /3i and {32 differ slightly between the two models. The skewed log-BS model yields the 
highest value of the log-likelihood function and smallest values of the AIC, BIC and HQIC statistics. 
From the values of these statistics, the skewed log-BS model outperforms the BS model and should 
be prefered. The likelihood ratio (LR) statistic to the null hypothesis A = is in accordance with the 
information criteria (LR = 5.88 and the associated critical level of the Xi at 5% is 3.84). 

In what follows, we shall apply the generalized leverage and local influence methods developed 
in the previous sections for the purpose of identifying influential observations in the skewed log-BS 
regression model fitted to the data set. Figure [21 gives the \d^iJ\ corresponding to 6 for different 
perturbation schemes and the generalized leverage. An inspection of Figure [21 reveals that based on 
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case-weight perturbation (Figure[2a)), we observed that the cases #2, #3, #10, #18 and #40 have more 
pronounced influence than the other observations. The case #1 appears with outstanding influence 
based on response perturbation (Figures Ob)). From the Figure [2tc) (covariate perturbation), the case 
#1 and #2 have more pronounced influence than the other observations. Figure [2td) reveals that the 
cases #9, #10 and #21 have influence on their own-fitted values. 

Based on Figure [2l we eliminated those most influential observations and refitted the skewed log- 
BS regression model. In Table[2]we have the relative changes of each parameter estimate, defined by 
RC = \{6j — and the corresponding SE, where 6'j(j) denotes the maximum likelihood esti- 

mate of 9j, after removing the ith observation. As can be seen, except for the case #21 corresponding 
to the parameter A, the relative changes for the maximum likelihood estimates of (32, a and A are very 
little pronounced. Also, the significance of these parameters are not modified in all cases considered. 
Case #21 represents the smallest value of the time to failure. Further, (3i becomes not significant in 
all cases considered similar to the skewed log-BS regression model fitted considering all observations 
(Table 

Table 2: Relative changes dropping the cases indicated, and the corresponding asymptotic standard 

errors. 

Pi (32 a A 



Dropping 


RC 


SE 


RC 


SE 


RC 


SE 


RC 


SE 


#1 


0.201 





180 


0.029 


1.614 


0.026 


0.341 


0.035 


0.543 


#2 


0.161 





180 


0.024 


1.618 


0.013 


0.345 


0.030 


0.544 


#3 


0.145 





180 


0.022 


1.619 


0.009 


0.347 


0.030 


0.544 


#9 


0.304 





181 


0.033 


1.665 


0.027 


0.369 


0.068 


0.609 


#10 


0.543 





180 


0.051 


1.675 


0.005 


0.358 


0.044 


0.588 


#18 


0.336 





177 


0.011 


1.558 


0.015 


0.347 


0.000 


0.569 


#21 


0.549 





132 


0.173 


1.125 


0.913 


0.150 


3.312 


0.406 


#40 


0.121 





176 


0.027 


1.626 


0.014 


0.346 


0.019 


0.565 



7 Concluding remarks 



In this paper we have introduced a log-Birnbaum-Saunders regression model with asymmetric errors, 
extending the usual log-BS regression model. The ran dom errors of the regression model follow 



LeivaetaL 



( 20101) . The estimation of the 



a skewed sinh-normal distribution, recently derived by 
model parameters is approached by the method of maximum likelihood and the observed information 
matrix is derived. We also consider diagnostic techniques that can be employed to identify influential 
observations. Appropriate matrices for assessing local influence on the parameter estimates under 
different perturbation schemes are obtained. The expressions derived are simple, compact and can be 
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easily implemented into any mathematical or s tatistical/econometric programmin g environr nent with 
nume rical linear algebra facilities, such as R (|R Development Core TeamL l2009h and Ox (jDoomik . 
20061) . among others, i.e. our formulas related with this class of regression model are manageable, 
and with the use of modern computer resources, may turn into adequate tools comprising the arsenal 
of applied statisticians. Finally, an application to a real data set is presented to illustrate the usefulness 
of the proposed model. 

As future research, it should be noticed that some generalizations of the proposed model could 
be done. For example, a skewe d log-BS regressio n model that allows us consider censored samples 
could be introduced. Following IXi and Weil (|2007h . one could introduce a skewed log-BS regression 
model in which the parameter a is considered different for each observation, i.e. to propose an het- 
eroscedastic skewed log-BS regression model. Also, a skewed log-BS nonlinear regression model 
could be proposed, and so forth. 
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Appendix 

After extensive algebraic manipulations, the quantities necessary to obtain the observed information 
matrix for the parameter vector Q = {/3^ , a, A)^ presented in the Section|2]are given by 



vAe) 



1 f 2 4 
4 2^- + ^ 



i2 



^K^2) 



hi = hi{0) 



6ie 



i2 



+ 



I £i2 

ACa / ii2<P{Ki2) >^^^fl^i2<P{Ki2 



a 



4 [ 

2 £-2 



4 t $(Ae.2) 

A f U<P{K^2) \%a(^2(t^{\U) 

2a\ $(Aei2) $(A62) 



^(^lCt>{K^2f 
^{H^2f 
A6l^»20(A^^2) 
^{Ki2? 



+ 



ACA f 620(A^.2) >?^U^2(t^{\i^2) A^fi0(A^,2) 



4 i $(Ae 



i2 



i2 



l[ 6l0(A62) >^%lil4>{K^2) 

2\ $(Ae.2) $(A62) 



HH^2) 

Ae.i620(Ae 



i2) 



<^{K^2y 



12 



ka(0) = A - ^ 



y |6i62 - ^1 + 2$(A62) + ^ 



g 

2 



i2 



+ o (Cjl + ^12) 



i2 



2c^ 

a^efj a$(A62) 
Aei20(Ae.2) 

2$(Aei2) la' 2 J ' a2$(A62) 
A^Ca^a0(A^j2) [ Ca^n \f ^c ^{^12) 

2$(Ae.2) I « ^ 2/1 ^^'+$(A62) 



_^ Aca0(A^j2) / in ^ C<xii2 ^ 



+ 



a^[\ii2) I a 



A^. 



i2 



0(Ae 



•i2; 



*(Ae 



i2j 



_Ci2 Cpi^ji 



; 7 ia\ e ^^2! , ACaA6i0(ACi2) 

Ca J 2^ji^j2 , Ca,^2 I t2\ 2c, 



0(A62 



Aca0(A^^2) / 61 Ca^j2 

2$(Ae.2) la 2 J $(A62 



_Ci2 CgCa 



^^C\ii\4>{Ki2) \ ii2 '/'(A62) 

2$(A62) I a ^ 2 /l^^'^$(A62) 

A620(A(,2) / e.2 CaCaU.^ , 0(Ae^\ 

«$(Aei2) I a ^ 2 i r^^'^$(Aei2)J' 



A;i3(0) 



(. (. ii2 

2 1^^^^^^-^;: 



+ 



$(Ae.2) ' 2$(Ae^2) 
if ..2 ^ A _ 1 ^ ^11 ^ Acie.20(A62) 
4r^" + a^ ^e^J^ 4$(A62) 

A^CA6lCi20(A62) , Aca6i\ ^C\ii\'^{Ki2f \ c , -^Ca^I 

Aa0(A6: 



e.2 + 



2 J ^<±'l^/V(;j2 
2 



2$(Ae.2^2 



Aca6i 1 620(ACi2)' 



^^2 + 



Aca6 



2$(Ae.2) r 2 J 2$(A62)2 ■ 2 

, n. Also, 

/oo 
w^(4 + Q;2w^)"^/VM$(Aw)(iw, 
-00 

/oo 
smh~^(Q;i(;/2)0(i(;)0(Ai(;)di(;, 
-00 

/oo 
-00 



13 



References 



Azzaline, A. (1985). A class of distributions which includes the normal ones. Scandinavian Journal of Statistics 

12, 171-178. 

Balakrishnan, N., Leiva, V., Lopez, J. (2007). Acceptance sampling plans from truncated life tests from gen- 
eralized Birnbaum-Saunders distribution. Communications in Statistics - Simulation and Computation 36, 
643-656. 

Bhatti, C.R. (2010). The Birnbaum-Saunders autoregressive conditional duration model. Mathematics and 
Computers in Simulation 80, 2062-2078. 

Bimbaum, Z.W., Saunders, S.C. (1969a). A new family of Ufe distributions. Journal of Applied Probability 6, 
319-327. 

Birnbaum, Z.W., Saunders, S.C. (1969b). Estimation for a family of life distributions with applications to 
fatigue. Journal of Applied Probability 6, 328-377. 

Cancho, V.G., Ortega, E.E.M., Paula, G.A. (2010). On estimation and influence diagnostics for log-Bimbaum- 
Saunders Student-t regression models: Full Bayesian analysis. Journal of Statistical Planning and Inference 
140, 2486-2496. 

Chan, P.S., Ng, H.K.T., Balakrishnan, N., Zhou, Q. (2008). Point and interval estimation for extreme-value 
regression model under Type-II censoring. Computational Statistics ans Data Analysis 52, 4040-4058. 

Castillo, N.O., Gomez, H.W., Bolfarine, H. (2009). Epsilon Birnbaum-Saunders distribution family: properties 
and inference. Statistical Papers. DOI : 10 . 1007/s00362-009-0293-x. 

Cook, R.D. (1986). Assessment of local influence (with discussion). Journal of the Royal Statistical Society B 
48, 133-169. 

Cox, D.R., Hinkley, D.V. (1974). Theoretical Statistics. London: Chapman and Hall. 

Desmond, A.F. (1985). Stochastic models of failure in random environments. Canadian Journal of Statistics 

13, 171-183. 

Desmond, A.F. (1986). On the relationship between two fatigue-Ufe models. IEEE Transactions on Reliability 
35, 167-169. 

Desmond, A.F., Rodriguez-Yam, G.A., Lu, X. (2008). Estimation of parameters for a Birnbaum-Saunders 
regression model with censored data. Journal of Statistical Computation and Simulation 78, 983-997. 

Diaz-Garcia, J. A., Leiva, V. (2005). A new family of life distributions based on the elhptically contoured 
distributions. Journal of Statistical Planning and Inference 128, 445^57. 

Doomik, J.A. (2006). An Object-Oriented Matrix Language - Ox 4, 5th ed. Timberlake Consultants Press, 
London. 



14 



Galea, M., Leiva, V., Paula, G.A. (2004). Influence diagnostics in log-Bimbaum-Saunders regression models. 
Journal of Applied Statistics 31, 1049-1064. 

Gomes, H.W., OUvares-Pacheco, J.F., Bolfarine, H. (2009). An extension of the generalized Bimbaum- 
Saunders distribution. Statistics and Probability Letters 79, 331-338. 

Guiraud, P., Leiva, V., Fierro, R. (2009). A non-central version of the Bimbaum-Saunders distribution for 
reliability analysis. IEEE Transactions on Reliability 58, 152-160. 

Kundu, D., Kannan, N., Balakrishnan, N. (2008). On the function of Birnbaum-Saunders distribution and 
associated inference. Computational Statistics and Data Analysis 52, 2692-2702. 

Leiva, V., Barros, M.K., Paula, G.A., Galea, M. (2007). Influence diagnostics in log-Bimbaum-Saunders re- 
gression models with censored data. Computational Statistics and Data Analysis, 51, 5694—5707. 

Leiva, V., Barros, M., Paula, G.A., Sanhueza, A. (2008). GeneraUzed Bimbaum-Saunders distributions applied 
to air pollutant concentration. Environmetrics 19, 235-249 . 

Leiva, V., Sanhueza, A., Angulo, J.M. (2009). A length-biased version of the Birnbaum-Saunders distribution 
with application in water quality. Stochastic Environmental Research and Risk Assessment 23, 299-307. 

Leiva, V., Vilca, P., Balakrishnan, N., Sanhueza, A. (2010). A skewed sinh-normal distribution and its properties 
and application to air pollution. Communications in Statistics - Theory and Methods 39, 426-443. 

Lemonte, A.J., Cordeiro, G.M. (2009). Bimbaum-Saunders nonlinear regression models. Computational 
Statistics and Data Analysis 53, 4441-4452. 

Lemonte, A.J., Cribari-Neto, F., Vasconcellos, K.L.P. (2007). Improved statistical inference for the two- 
parameter Bimbaum-Saunders distribution. Computational Statistics and Data Analysis 51, 4656^681. 

Lemonte, A.J., Ferrari, S.L.P., Cribari-Neto, F. (2010). Improved likelihood inference in Bimbaum-Saunders 
regressions. Computational Statistics and Data Analysis 54, 1307-1316. 

Lemonte, A.J., Simas, A.B., Cribari-Neto, F. (2008). Bootstrap-based improved estimators for the two- 
parameter Bimbaum-Saunders distribution. Journal of Statistical Computation and Simulation 78, 37^9. 

Lepadatu, D., Kobi, A., Hambli, R. and Barreau, A. (2005). Lifetime multiple response optimization of metal 
extrusion die. Proceedings of the Annual Reliability and Maintainability Symposium, 31-42. 

Marshall, A.W., Olkin, I. (2007). Life Distributions. Springer: New York. 

McCool, J.I. (1980). Confidence limits for WeibuU regression with censored data. IEEE Transactions on Reli- 
ability 29, 145-150. 

Nocedal, J., Wright, S.J. (1999). Numerical Optimization. Springer: New York. 

Press, W.H., Teulosky, S.A., Vetterling, W.T., Flannery, B.P. (2007). Numerical Recipes in C: The Art of Scien- 
tific Computing, 3rd ed. Cambridge University Press. 



15 



R Development Core Team (2009). R: A Language and Environment for Statistical Computing. Vienna, Austria. 

Rieck, J.R. (1989). Statistical Analysis for the Bimbaum-Saunders Fatigue Life Distribution. Ph.D. disserta- 
tion, Clemson University. 

Rieck, J.R., Nedelman, J.R. (1991). A log-linear model for the Bimbaum-Saunders distribution. Technometrics 
33,51-60. 

Saunders, S.C. (1974). A family of random variables closed under reciprocation. Journal of the American 
Statistical Association 69, 533-539. 

Tisionas, E.G. (2001). Bayesian inference in Bimbaum-Saunders regression. Communications in Statistics - 
Theory and Methods 30, 179-193. 

Xi, EC, Wei, B.C. (2007). Diagnostics analysis for log-Bimbaum-Saunders regression models. Computational 
Statistics and Data Analysis 51, 4692^706. 

Xiao, Q., Liu, Z., Balakrishnan, N., Lu, X. (2010). Estimation of the Bimbaum-Saunders regression model 
with current status data. Computational Statistics and Data Analysis 54, 326-332. 

Xu, A., Tang, Y. (2010). Reference analysis for Bimbaum-Saunders distribution. Computational Statistics and 
Data Analysis 5A, 185-192. 

Wei, B.C., Hu, Y.Q., Fung, W.K. (1998). Generahzed leverage and its applications. Scandinavian Journal of 
Statistics 25, 25-37. 

Wu, J., Wong, A.C.M. (2004). Improved interval estimation for the two-parameter Bimbaum-Saunders distri- 
bution. Computational Statistics and Data Analysis 47, 809-821. 



16 



